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We use discontinuous molecular dynamics and grand-canonical transition-matrix Monte Carlo 
simulations to explore how confinement between parallel hard walls modifies the relationships be- 
tween packing fraction, self-diffusivity, partial molar excess entropy, and total excess entropy for 
binary hard-sphere mixtures. To accomplish this, we introduce an efficient algorithm to calculate 
partial molar excess entropies from the transition-matrix Monte Carlo simulation data. We find 
that the species-dependent self-diffusivities of confined fluids are very similar to those of the bulk 
mixture if compared at the same, appropriately defined, packing fraction up to intermediate values, 
but then deviate negatively from the bulk behavior at higher packing fractions. On the other hand, 
the relationships between self-diffusivity and partial molar excess entropy (or total excess entropy) 
observed in the bulk fluid are preserved under confinement even at relatively high packing fractions 
and for different mixture compositions. This suggests that the partial molar excess entropy, calcu- 
lable from classical density functional theories of inhomogeneous fluids, can be used to predict some 
of the nontrivial dynamical behaviors of fluid mixtures in confined environments. 



I. INTRODUCTION 

Confinement of a fluid can substantially modify its 
physical properties. 1-3 Although the associated changes 
to thermodynamics and structure can often be predicted 
by theory for simple fluids, estimating the implications 
of confinement for transport coefficients is more challeng- 
ing. Since the latter remains key to understanding many 
systems of scientific and technological interest, the dis- 
covery of even simple heuristics could have significant 
impact. 

In this spirit, one useful line of inquiry is to use molecu- 
lar simulations to understand the behavior of simple fluid 
models in controlled confinement conditions. Here, we 
study hard-sphere (HS) fluid mixtures, both in bulk and 
as thin films confined between smooth hard walls. The 
HS model captures much of the important fundamental 
physics of fluids, most notably entropic packing effects 
associated with excluded volume interactions. While it 
is arguably the most basic fluid model, its historically sig- 
nificant role in understanding the bulk fluid-solid phase 
transition, 4,5 its structural similarity to simple atomic 
liquids, 6 and its ability to quantitatively describe exper- 
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iments on colloidal suspensions, 7 9 make it the standard 
reference system for fluids. 

Given the simplicity of the HS potential, it is per- 
haps surprising that numerically precise and compre- 
hensive data for the monatomic HS fluid confined be- 
tween parallel hard walls have been obtained only very 
recently. 10 These data show that, for a wide range of 
confinement conditions, the fluid behavior is similar to 
that of the bulk, provided that the comparison is made 
using an appropriately defined density or packing frac- 
tion. Specifically, the density should be defined based on 
the total volume of the system as opposed to the smaller 
volume accessible to fluid particle centers. 10,11 The ther- 
modynamic and kinetic behavior of the HS fluid con- 
fined between hard walls starts to deviate significantly 
from the bulk behavior only at relatively high particle 
densities and restrictive pore sizes. 10 Based on compar- 
ison between the confinement shifted solid-fluid phase 
boundaries as calculated by Fortini et al., 12 and devi- 
ations from bulk behavior in quantities like transverse 
diffusivity, pressure, etc. it appears that the origin of 
these effects is packing frustration due to the confined 
geometry being incommensurate with the formation of 
an integer number of fluid layers between the hard walls. 

A second observation based on recent simulations of 
the monatomic HS fluid is that the relationship be- 
tween excess entropy (relative to an ideal gas) and self- 
diffusivity is nearly identical for bulk and confined sys- 
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terns. 10,11 In fact, excess entropy is a significantly more 
accurate predictor than average density for how confine- 
ment of the HS fluid affects its self-diffusivity. A similar 
observation also holds for predicting the effects of con- 
finement for the dynamics of the equilibrium Lennard- 
Jones and square- well fluids. 13 The use of excess entropy 
to estimate the self-diffusivities of confined fluids was 
originally motivated by earlier simulations which demon- 
strated a strong correlation between the two quantities 
for a number of bulk fluids. 14-20 More recent simulations 
indicate that excess entropy (and its two-body approxi- 
mation based on the pair correlation function 2 1-23 ) can 
also provide new insights into the connections between 
the static structure, thermodynamics, and dynamics of 
supercooled liquids. 24-27 

Can the relationships for monodisperse fluid systems 
discussed above be generalized for multicomponent sys- 
tems? From a practical perspective, understanding the 
behaviors of equilibrium mixtures will have benefits, 
since many experimental systems are inherently poly- 
disperse. There is also an urgent need for improved 
predictive capabilities for transport coefficients in fluid 
mixtures due to their relevance in industrial separation 
processes. 28 Finally, understanding equilibrium mixtures 
will serve as a foundation for future studies of confined, 
supercooled liquids. There are a number of fundamen- 
tal open questions regarding these latter systems, in- 
cluding predicting how confinement will shift the glass 
transition. 29-31 The reason that mixtures are important 
for studying supercooled, confined liquids is that their 
polydispersity frustrates the crystallization that other- 
wise rapidly occurs in monodisperse samples. 

In this paper, we use molecular dynamics and 
grand-canonical transition-matrix Monte Carlo 
simulations 32 ' 33 to explore whether density, excess 
entropy, or partial molar excess entropy can be used 
to forecast how confinement will affect the species- 
dependent self-diffusivity of hard-sphere mixtures. One 
methodological outcome of our investigation is the 
introduction of an efficient algorithm to calculate partial 
molar excess entropies of mixtures from grand-canonical 
transition-matrix Monte Carlo simulation data. Our 
main scientific finding is that previous results concerning 
the relationship of excess entropy and self-diffusivity 
of confined monodisperse systems appear to generalize 
well to the fluid mixtures we investigate here. This 
suggests that knowledge of the behavior of bulk mixtures 
together with predictions for how confinement modifies 
their thermodynamics (from, e.g., density functional the- 
ory) can be used to accurately predict the single-particle 
dynamics of confined mixtures. 



II. SIMULATION AND THEORETICAL 
METHODS 

In this section, we describe the binary HS fluid in- 
vestigated here and methods used to calculate its exact 



thermodynamic and kinetic properties. 



A. Model and Simulations 

We studied a binary mixture of HS particles with 
disparate particle diameters and masses. The ratio of 
the particle diameters is given by a\jai = 1.3 and the 
particle masses are proportional to their volume, i.e. 
toi/to2 = <j\I<J%- We have chosen these parameter values 
to mimic recent experiments on binary colloidal mixtures 
in confinement. 34 This choice sets the stage for future 
studies in which we plan to make comparisons to exper- 
imental observations in the supercooled fluid. 

We calculated the thermodynamic properties of the 
bulk and confined HS mixture using grand-canonical 
transition-matrix Monte Carlo (GC-TMMC). 35 ' 36 These 
simulations are conceptually equivalent to a series of 
semigrand simulations performed over a range of fluid 
densities stitched together using ghost insertion/deletion 
moves. The reader is referred to earlier work 35,36 for more 
details. The primary quantity yielded by GC-TMMC 
is the particle number probability distribution. For a 
multicomponent system, this distribution is multidimen- 
sional and is denoted as n(7V; fi, V, T). Mathematically, 
n(iV; /x, V, T) represents the probability of observing the 
particle number vector N = (Aq, A 2 , ...Ni, ...), where N t 
is the number of particles of species i, in a system of vol- 
ume V at temperature T and imposed set of chemical 
potentials /x = /12, ■■■Hi, •••), where ^ is the chemical 
potential of species i. 

We simulated the bulk and confined binary HS mix- 
ture with specified activities £1 = £2 = 7.8446 and 
(3 = 1.0, where & = exp(/3^)/A 3 , (3 = l/(k B T) {k B 
is the Boltzmann constant), T is temperature, and Aj 
is the de Broglie wavelength of species i. Wc required 
each (Aq, A2) pair to be visited by the system a mini- 
mum of two million times. We used a system volume of 
V = I25er 3 for the bulk fluid simulations, and we em- 
ployed smooth hard walls of surface area of 49er 2 for the 
confined fluid simulations. In addition, to verify that 
system-size effects were negligible, we performed a series 
of simulations employing at least twice the volume. For 
the slit-pore geometry, we calculated density profiles at 
each (Ai, A 2 ) pair using one hundred bins. 

We used discontinuous molecular dynamics (DMD) 
simulations 37 to calculate the species-specific self- 
diffusivity of the binary system. For the bulk fluid, we 
employed a cubic simulation cell of V, and we applied 
periodic boundary conditions in all three directions. For 
the confined fluid, we used a rectangular parallelepiped 
simulation cell of V = H x H y H z , where H z is the to- 
tal (i.e., not the center-accessible) distance between the 
confining surfaces. We applied periodic boundary condi- 
tions in the x and y directions and we placed perfectly 
reflecting, smooth hard walls so that particle surfaces 
were trapped in the region < z < H z . We obtained 
the self-diffusivity Di of the individual fluid components 
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by fitting the long-time (t ^> 1) behavior of the aver- 
age mean-squared displacement of the zth particles to 
the Einstein relation (Ar|). = 2dDit, where Ar^ corre- 
sponds to the mean-square displacement per particle for 
i type particles (i = 1 or 2) in the d periodic directions 
(d = 2,3 for the confined and bulk fluid, respectively). 
We used TV = 2182 particles for all DMD simulations. 
We also ran simulations with TV = 1091 and 4364 for 
several state points, and we found finite-size effects to be 
within statistical uncertainty of the simulation data. 



B. Total Excess Entropy S e 



(4) provides the total entropy of the system in terms of 
information obtained from GC-TMMC simulations and 
is completely general. Since this work deals with HS 
systems, in what follows, the discussion is restricted to 
mixtures of spherical particles. 

We study both homogeneous (bulk) and inhomoge- 
neous fluid mixtures in this work. We first consider the 
ideal gas reference state of a homogeneous fluid. The 
Helmholtz free energy F lg of a bulk ideal gas mixture is 



f3F ig {N,V 7 T) = -J2 ln 



V 



3Ni 



(5) 



Knowledge of the particle number distribution 
U(N;fj,,V,T), which can be obtained from GC-TMMC 
simulation, is tantamount to knowledge of the system's 
free energy as a function of density and composition. One 
can show that the Helmholtz free energy F of a system 
is related to the particle number distribution 



(3F(N, V, T) 



ln n(iV;/x,y,T) 

n n(0;M,V,T) 



(3j2^N t (1) 



where n(TV; p, V, T) is the normalized particle num- 
ber distribution obtained from a GC-TMMC simulation. 
While IT(TV; p, V, T) is the primary quantity of interest, 
additional quantities can also be determined during the 
course of a GC-TMMC simulation. 38 For example, in this 
work, for a given particle number vector TV, we calcu- 
lated the average total potential energy [/(TV), and, in 
the case of the slit pore geometry, the average density 
profile of each species pi(r;N). From thermodynamics, 
it is known that 



F = E-TS (2) 



where S is the entropy and E is the total energy 



E = K + U (3) 



where K is the kinetic energy and U is the potential 
energy. It is straightforward to show that 



S(N,V,T) , Il(N;p,V,T) 



In- 



n(0;M,V,T) 



+ pE(N)-(3j2»i N i 



(4) 

To calculate the partial molar excess entropy s° x , it is 
first necessary to calculate the excess entropy 5 CX . In this 
work, S cyi is the entropy difference between the system 
of interest and an ideal gas with the same temperature 
and spatial distribution of density and composition. Eq. 



Substituting Eqs. (2) and (3) into Eq. (5), an expres- 
sion for the ideal gas entropy can be obtained. It then 
immediately follows that the excess entropy is 



S° X (JV, V, T) 
k B 



ln n(Tv ; ^y,T) 



I3U(N) 



n(0;€,V,T) 
^-TViln^+lnTVJ-TV^lnTVi 

i 

+TV i lnp j 



(6) 



where pi is the bulk number density of species i. 

Now we consider an inhomogcneous fluid mixture. In 
this case, the free energy of the system also depends on 
the density profile of each species in the system. In other 
words, the free energy is a functional of the set of den- 
sity profiles p — {pi(r) 7 p 2 (r), ...pi(r), ...}. 39 The excess 
entropy can be calculated in a manner similar to the ho- 
mogeneous case, where the total entropy of the system S 
is again given by Eq. (4). The ideal gas Helmholtz free 
energy of an inhomogeneous fluid mixture is 40,41 

pF i °([p],N,V,T)= f drJ2^(r){^[ Pt (r)Af] - 1} . 
J i 

(7) 

Making use of Eqs. (2), (3), (4), and (7), one can show 
that the excess entropy of an inhomogeneous fluid mix- 
ture is given by an expression similar to that of the ho- 
mogeneous fluid 



S-jNV,T) = ln fL(N;S,V,T) 
k B n(0;€,V,T) 

~Ni ln& + In TVJ - TV, In TV; 

i 

+ J dr Pi (r;N)ln Pi (r;N) (8) 



where pi(r;N) is the density profile of species i for a 
specified particle number pair (TVi,TV 2 ). Eq. (8) is the 
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multicomponcnt extension of the expression derived ear- 
lier by Mittal et a/. 11 In addition, notice that when Pi(r) 
is uniform, the expression for the total excess entropy of 
a homogeneous fluid mixture Eq. (6) is recovered. 

C. Partial Molar Excess Entropy sf x 

Given the excess entropy S cx of a mixture, the partial 
molar excess entropy s® x of component i is defined as 42 

V WV « / T,p,N m 

Notice that the derivative is taken at fixed temperature 
T, pressure p, and number of each species j other than 
i. However, the expressions for calculating the excess 
entropy using a GC-TMMC simulation, Eqs. (6) and (8), 
are functions of volume, not pressure. Therefore, it is not 
straightforward to take the partial derivative with respect 
to Ni directly while fixing the pressure. To circumvent 
this difficulty, we first express the total differential of S cx 
in terms of T, p, and N 

Imposing the constraint of fixed T, V, and Njui, the par- 
tial molar excess entropy can be alternatively expressed 
as 

f dS cx \ ( ds cx \ ( dp \ 

sT = {-m) TyNj[i] - N \~dp-) T , x {m) T , VtNm 

(11) 

where N = J2iN, s cx = S cx /N is the excess entropy 
per particle, and x = (x\,X2, ■■■Xi, ...) represents the mole 
fraction vector for the species. 

Eq. (11) is a general expression and provides the 
framework for calculating the partial molar excess en- 
tropy using information obtained from a GC-TMMC sim- 
ulation. In particular, the partial derivative taken at 
fixed composition suggests that a natural way to calcu- 
late the partial molar excess entropy is along an isoplcth 
while varying the fluid density. 

To evaluate the partial derivatives in Eq. (11), we 
adopt an approach where we construct numerically the 
required functions (e.g., s cx (p) at fixed T and x) and 
then fit them to polynomials, whose derivatives can be 
evaluated analytically. However, from a numerical per- 
spective, this is awkward to implement directly from raw 
simulation data because N can only take on integer val- 
ues, and thus any function of Ni is discontinuous. One 



possible route to circumvent this difficulty involves work- 
ing in terms of the mean or ensemble-averaged analogs 
of the quantities in Eq. (11), since mean values can take 
on a continuous range of values. In particular, instead of 
Ni, we use the mean particle number (Ni), which is 

(N i ) = , £ i N i n(N;t,V,T). (12) 

N 

It is likewise more convenient to work with the average 
excess entropy (S ex ), 

(S cx ) = J2 S cx (iV, V, T) fl(N; £, V, T) (13) 

JV 

where S CX (N ,V,T) is given by Eq. (6) or (8). Finally, 
the pressure p, which does not require any averaging, is 
simply 

ppV=-lnfl(0;£,V,T) . (14) 

Notice that the mean quantities and the pressure are ex- 
plicit functions of the activities or chemical potentials. 
Given the particle number probability distribution ob- 
tained from TMMC simulation at specified the dis- 
tribution can be determined at other activities by using 
histogram reweighting, 43 thus allowing for the calcula- 
tion of (Ni), (S cx ), and p as a function of £, at constant 
volume and temperature. Our numerical strategy, which 
we describe below, is to construct the required functions 
indicated in Eq. (11) by varying £ and then fitting them 
to polynomials. 

We first focus on the calculation of the partial deriva- 
tive of the excess entropy per particle with respect to 
pressure at fixed composition and temperature. For 
the binary HS mixture, an initial activity of one of the 
species, say species 1, is set to some arbitrary value, usu- 
ally corresponding to a dilute vapor. The activity of the 
other species is then adjusted such that the average com- 
position of the fluid corresponds to the desired isopleth 
(within a fractional tolerance of 1 x 1CP 8 ), and the av- 
erage excess entropy per particle and pressure for this 
activity pair are calculated and recorded. The activity of 
species 1 is then incremented and the process is repeated 
until total densities ranging from vapor-like to liquid-like 
values are explored. This iterative process yields the av- 
erage excess entropy per particle as a function of pressure 
along an isopleth. In the bottom panel of Figure 1, we 
plot (s cx ) versus p for the bulk fluid mixture along an 
isoplcth corresponding to equal volume fraction of each 
species. Notice that the curve is smooth and continuous, 
allowing for the straightforward numerical determination 
of its derivative. 

We now focus on the calculation of the terms in Eq. 
(11) involving partial derivatives at fixed average particle 
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FIG. 1: Examples of quantities needed to calculate the par- 
tial molar excess entropy. In the top and middle panel, we 
show the average excess entropy per particle and pressure, 
respectively, as a function of the average number density 
((/**) = (Ni) /V) of species 1 at fixed average number den- 
sity of species 2. The average excess entropy per particle as 
a function of pressure is shown in the bottom panel at a fluid 
composition corresponding to equal volume fractions of each 
species. 



number. Conceptually, this is done by numerically con- 
structing (S ex ) and p as a function of (iVj) while holding 
the other average particle numbers fixed. The procedure 
is similar to that used to construct the (s ex )-p curve. 
Consider, for illustrative purposes, the situation where 
(S ex ) and p are sought as a function of (Nx) at fixed av- 
erage number of species 2, say (N%). To do this, an ini- 
tial activity £i is specified, and the value of £2 is adjusted 
such that (N2) = (N%) (within a fractional tolerance of 
1 x 10 -8 ). It should be noted that simply fixing £2 does 
not necessarily fix (N2). The quantities (S ex ), p, and 
(Nx) are then calculated and recorded for this pair of ac- 
tivities, and a new value of £1 is specified. This process is 
repeated until the required quantities are obtained over 
the desired range of (Nx) values. In the top and middle 
panels of Figure 1, examples of these curves are shown. 
Again, notice that the data are smooth and continuous. 
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FIG. 2: Partial molar excess entropy and average excess 
entropy per particle as a function of total fluid density for the 
bulk hard-sphere mixture at a composition corresponding to 
equal volume fractions of each species. 

From a practical point of view, one only needs to con- 
struct the portions of these curves that coincide with the 
isopleth used to determine the (s ex )-p curve. 

Using the procedure outlined above, the partial molar 
excess entropy sf x can be calculated as a function of total 
density at fixed composition. In Figure 2, we plot the 
partial excess entropy of both species as a function of 
average total density (p) = (N) /V for the bulk fluid at 
equal volume fraction of each species. Also plotted in the 
same figure is the average excess entropy per particle. As 
a test of thermodynamic self-consistency, we have also 
verified that the average excess entropy per particle is 
equal to the weighted sum of the partial molar excess 
entropies. 



III. RESULTS AND DISCUSSION 

A. Relationship between packing fraction and 
self-diffusivity 

First, we compare the self-diffusivity Di of bulk and 
confined fluid mixtures as a function of packing fraction 
(j). In this subsection, the goal is to see if packing fraction 
alone can describe the individual component diffusivity 
for a confined fluid mixture. Before presenting the sim- 
ulation data, we want to reiterate here our earlier pro- 
posal in recent studies that the density of confined fluids 
should be defined based on the total volume as opposed 
to the particle center accessible volume. 10,11 ' 13 When cal- 
culated in this manner, we found that packing fraction 
described the self-diffusivity of a number of pure simple 
fluids in bulk and confinement. We refer the readers to 
Refs. 10,11,13 for further details. 



6 




FIG. 3: Self-diffusivity versus total packing fraction for binary 
HS fluid in bulk and confined between hard walls. 



FIG. 4: Self-diffusivity versus partial molar excess entropy 
for bulk and confined binary hard-sphere fluid between hard 
walls, (a) 01 = cj> 2 (b) 01 = 2. f 9702. 



The predictions for the self-diffusivity Di (i = 1, 2) 
versus total packing fraction <p for the binary HS mix- 
ture are presented in Fig. 3. Here, the total packing 
fraction is given by <j> — (71-/6) (pi a\ + Pio^). Data 
are presented for two mixture compositions, equal vol- 
ume fraction (0i = $2), and equal particle numbers 
{4>i = 2.19702), i-e., equal mole fraction. The top and 
bottom plots correspond to the diffusivity Di of particles 
i = 1 or 2 respectively and the symbols represent the con- 
finement data as shown in the legend in Fig. 3. It is clear 
from this plot that the diffusivities under confinement are 
very close to their bulk values up to intermediate pack- 
ing fractions (cf> < 0.4) but then deviate quantitatively at 
higher values of <j). The diffusivities of the confined fluid 
can be as small as 30% of their bulk values at the high- 
est packing fraction and smallest pore-size investigated 
in this work (0 = 0.45, and H — ia\). Stated differently, 
if one uses packing fraction as a predictive tool for the 
diffusivity at the high packing fractions, then one will 
greatly overestimate the actual value. It is still remark- 
able that for a wide range of conditions, the data for 
confined fluids fall very close to the bulk curve. 



B. Relationship between partial molar excess 
entropy and self-diffusivity 

Now, we investigate if a different thermodynamic quan- 
tity, i.e. partial molar excess entropy s| x , correlates with 
the confined fluid's self-diffusivity more reliably than 
packing fraction. Figure 4 shows the self-diffusivity data 
for both bulk and confined fluid mixtures plotted ver- 
sus sf x for the same conditions as Fig. 3. Note that the 
confined fluid data for a components 1 and 2 approxi- 
mately collapse onto the bulk curve over the whole range. 
This has potentially powerful implications for predict- 
ing the diffusivity of confined mixtures. Specifically, one 
can use existing thermodynamic theories, such as den- 
sity functional theory, for calculating the partial molar 
excess entropies for these mixtures in confinement and 
then use the bulk D — s? x relationships to predict dif- 
fusivities under confinement. The fact that confinement 
does not significantly change the relationships between 
self-diffusivity and partial molar excess entropy does not 
mean that these two quantities remain constant when 
the fluid is confined. Rather, it signifies that confine- 
ment affects these quantities in a way that preserves the 
relationship observed between the two for the bulk mix- 
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FIG. 5: Density profiles for components 1 and 2 for a binary 
hard-sphere fluid confined between hard walls. 



ture. In fact, we do observe confinement induced effects 
in our simulations. For example, changes in the self- 
diffusivity due to confinement can be seen in Fig. 3. This 
is also accompanied by noticeable changes in fluid struc- 
ture, specifically the formation of fluid layering, which 
can be seen in the density profiles presented in Fig. 5. 



C. Relationship between total excess entropy and 
self-diffusivity 

Here, we explore if the excess entropy itself can be used 
in combination with bulk fluid behavior to predict single- 
particle dynamics. Figure 6 shows the self-diffusivity 
data for (a) 4>\ — 4>2 and (b) 4>i — 2.197^2 in bulk (lines) 
and under confinement (symbols). For these very differ- 
ent compositions, one can see that the self-diffusivities 
collapse onto a species -specific curve independently of 
pore width. 

To see if the above findings hold across the entire com- 
position range of the mixture, we have also generated dif- 
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FIG. 6: Self-diffusivity versus total excess entropy for binary 
hard-sphere fluid in bulk and confined between hard walls. 
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FIG. 7: Self-diffusivity versus (a,b) total excess entropy, and 
(c,d) partial molar excess entropy for binary hard-sphere fluid 
in bulk and confined between hard walls. The particle mole 
fraction values are as shown in the legend. 



fusivity and excess entropy data at other compositions. 
These additional results are given in Fig. 7. In Figs. 7a 
and 7b, we find a remarkable collapse of data when the 
self-diffusivity is plotted against the excess entropy. For 
each species, we now find that the Di — s cx relationship 
holds, to a very good approximation, independently of 
pore width and composition. For completeness, the self- 
diffusivity is plotted against the partial excess entropy in 
Figs. 7c and 7d. In this case, the data collapse indepen- 
dently of pore width but not composition. 



IV. CONCLUSIONS 

In this paper, we present a systematic investiga- 
tion of the relationships between packing fraction, self- 
diffusivity, partial molar excess entropy, and total excess 
entropy for a binary HS mixture with components differ- 
ing in their sizes and masses. To a very good approxi- 
mation, the same relationship between the self-diffusivity 
and the packing fraction in both bulk and confined fluids 
is obeyed up to intermediate packing fractions. However, 
the deviations from this relationship at higher packing 
fractions and smaller pore sizes can be as high as 70%, 
indicating that the behavior of the confined fluid under 
such conditions differs significantly from the bulk. We 
find that the excess entropy and partial molar excess 
entropy are quantitatively more accurate predictors of 
single-particle dynamics (i.e., self-diffusivity) under con- 
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finement than total packing fraction. This conclusion is 
based upon simulation data generated over a broad range 
of compositions. 

Another outcome of this study is the introduction of 
a method to calculate the partial molar excess entropy 
from GC-TMMC simulations and histogram rewcighting. 
The method can be generalized to different fluid mixtures 
and fluid models. We are currently exploring the behav- 
ior of fluids with attractive interactions in the presence 
of attractive or repulsive surfaces as well as the behavior 
of the presently studied confined binary HS mixture in 
its supercooled state. In future studies, we also plan to 
use density functional theory to calculate the thermody- 
namics of mixtures, which will enable us to use the rela- 
tionships between self-diffusivity and excess entropy pre- 
sented here to make predictions about the single-particle 
dynamics of confined mixtures for a much broader range 
of conditions. 
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